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We study the f{R) theory of gravity using metric approach. In particular we investigate the 
recently proposed model by Hu-Sawicki, Appleby — Battye and Starobinsky. In this model, the 
cosmological constant is zero in flat space time. The model passes both the Solar system and the 
laboratory tests. But the model parameters need to be fine tuned to avoid the finite time singularity 
. . . recently pointed in the literature. We check the concordance of this model with the //(z) and baryon 

00 ' acoustic oscillation data. We find that the model resembles the ACDM at high redshift. However, for 

, some parameter values there are variations in the expansion history of the universe at low redshift. 
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It is remarkable that different data sets of complementary nature such as supernovae, baryon oscillations, galaxy 
clustering, microwave back ground and weak lensing all taken together strongly support the late time acceleration 
O i' of universe. In the standard lore, one assumes that the history of universe is described by the general relativity 
D , (GR). The late time acceleration can easily be captured in this frame work by introducing a scalar field with large 
, negative pressure known as dark energy In view of the fine tuning problem, the scalar field models, specially those 
with tracker like solutions, are more attractive compared to the models based on cosmological constant. At present, 
PsJ ■ observations are not in a position to reject or to establish the dark energy metamorphosis. A host of scalar field 
J> ' models have been investigated in the literature. The scalar field models can fit the data but lack the predictive power. 
, It then becomes important to seek the support of these models from a fundamental theory of high energy physics. 

■ One can question the standard lore on fundamental grounds. We know that gravity is modified at small distance 
scales; it is quite possible that it is modified at large scales too where it has never been confronted with observations 
directly. It is therefore perfectly legitimate to investigate the possibility of late time acceleration due to modification 
' of Einstein-Hilbert action. It is tempting to study the string curvature corrections to Einstein gravity amongst 
, which the Gauss-Bonnet correction enjoys special status. A large number of papers are devoted to the cosmological 
L| implications of string curvature corrected gravity^, 0, [1, @, 0, M, S [l^j Ef- These models suffer from several 
. ^ problems. Most of these models do not include tracker like solution and those which do are heavily constrained by the 
thermal history of universe. For instance, the Gauss-Bonnet theory with dynamical dilaton might cause transition 
5h from matter scaling regime to late time acceleration allowing to alleviate the fine tuning and coincidence problems. 
. . . 1 However, it is difficult to reconcile this model with nucleosynthesis Another possibility of large scale modification 
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is provided by non-locally corrected gravity which typically involves inverse of d'Alembertian of Ricci scalar. The 
non-local construct might mimic dark energy; the model poses technical difficulties and there has been a little progress 
in this direction [l2|. The large scale modification may also arise in extra dimensional theories like DGP model which 
contains self accelerating branc. Apart from the theoretical problems, this model is heavily constrained by observation. 



On purely phenomenological grounds, one could seek a modification of Einstein gravity b y re placing the Ricci scalar by 
/(i?). The f{R) gravity models have been extensively investigated in past five years [l3l [l4l. [la. fig . The /(i?) gravity 
theories giving rise to cosmological constant in low curvature regime are plagued with instabilities and on observational 
grounds they are not distinguished from cosmological constant. The recently introduced models of f(R) gravity by 
Hu-Sawicki and Starobinsky (referred as HSS models hereafter) with disappearing cosmological constant [13, EBl have 
given rise to new hopes for a viable cosmolo gica l model within the framework of modified gravity {f{R) gravity model 
with similar properties is proposed in Ref.[l^). These models contain Minkowski space time as a solution in the 
low curvature regime which is an unstable solution. In high curvature regime these models reduce to cosmological 
constant. Both the first and the second derivatives of f{R) with respect to R are positive. The positivity of the first 
derivative ensures that the scalar degree of freedom, a characteristic of any f(R) theory, is not tachyonic where as the 
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positivity of second derivative tells us that graviton is not ghost thereby guaranteing the stability. In Starobinsky 



parametrizationflSj, f{R) is given by, /(i?) = R + A 
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The HSS models can evade solar physics 



constraints provided that the model parameters are chosen properly. An important observation has recently been made 
by by Appleby — Battye and Forolov[2l], [S^l (see also [l^). The minimum of scalaron potential which corresponds 
to dark energy can be very near to = or equivalently R = oo. As pointed out in Ref.fS], the minimum should 
be near the origin for solar constraints to be evaded. Hence, it becomes most likely that we hit the singularity if the 
parameters are not fine tuned. 



In order to check whether the f{R) gravity theory is cosmological viable or not, it is necessary that this theory must 
be compatible with the observations. In this work, we study the Starobinsky model using the data from the recent 
observations which include H{z) , Hubble parameter at various red-shifts and the Baryon Acoustic Oscillation (BAO) 
peak from Sloan Digital Sky Survey (SDSS). 



This paper is organised as follows. In Section II, we describe the general properties of the model highlighting the 
fine tuning problem. The Friedmann equation and the special cases of trace equation is studied in Section HI. The 
cosmological constraints from the recent observations are described in Section IV. Finally, Section V contains the 
results and discussions. 



II. NATURALNESS OF THE MODEL 



In this section we shall revisit f(R) models with disappearing cosmological constant. These models have potential 
capability of being distinguished from ACDM and could lead to a viable cosmological model. However, even at the 
background level this class of models are fine tuned. In what follows we shall explicitly bring out these features 
specializing to Starobinsky parametrization. 



The action of f{R) gravity is given by[Tj|, 
which leads to the following equation of motion 



m_ 

IGttG 



-g d'^x, 



(1) 



(2) 



Here prime (') denotes the derivatives with respect to R. The f(R) gravity theories apart from a spin-2 object 
necessarily contain a scalar degree of freedom. Taking trace of Eq.® gives the evolution equation for the scalar 
degree of freedom. 



□/' = ^(2/-/'i?) + ^T. 
It would be convenient to define scalar function <j) as 

=./'-!, 

which is expressed through Ricci scalar once f(R) is specified. 



(3) 
(4) 



We can write the trace equation (equation ([3])) in the term of V and T as 



The potential can be evaluated using the following relation 

dV dV d(j) 1 



dR d6dR .^y-rmr. (6) 
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Recently Hu, Sawicki and Starobinsky proposed a functional form of /(i?) with the desirable properties of a viable 
model [I7|,ll8l. In this paper we shall consider the model in Starobinsky parametrization: 



f{R)=R + XRo 



R_ 



(7) 



Here n and A are greater than zero. And Rq is of the order of presently observed cosmological constant, A ~ SirGpuac- 
The properties of this model can be summarized as follows: 

1. In the absence of matter is always a solution, Muir^o f{R) = 0. However, as /" < 0, flat space time is 
unstable. 

2. For \R\ Rq, f{R) = R — 2A(oo). The high-curvature value of the effective cosmological constant is A(oo) = 
XRo/2. 

3. The stability conditions for the adopted model are 

fiR) > 0, f'{R) > 0. (8) 

In the Starobinsky model the scalar field (j) is given by 

2nXR 

-R) = 52 • (9) 

' i?o(l + 

We can compute V{R) for a given value of n. In case of n = 1, we have 

To " ^^^^^^{(-8-40?/-56/-24/)A+(32/ + lly3 + 2l2;5_3y^)A2} 
A^ 

— -tan~iy, (10) 



where y — R/ Rq. In case and of n = 2 we have, 
V 1 



{ (-160 - 1280y2 - 2880?;'* - 2560?;*^ - 800y*) A 



Ro ~ m{i + y^f 

+ (I05y + 595y^ + 215^ + 106?/^ + 595y^ + lOSy") A^} 
35A2 
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tan-iy. (11) 



In the FRW background, the trace equation (equation ([3])) can be rewritten in the convenient form 

•• •, dV SirG 

cj> + 3Hcl)+ — ^—p. (12) 

The time-time component of the equation of motion ([2]) gives the Hubble equation 

d(ln/')„ lf~f'R SttG 
H^ + ^ — U.H+-- ^ — = p. (13) 

We recover Einstein gravity in the limit /' = 1 . The simple picture of dynamics which appears here is the following: 
above infrared modification scale (Rq), the expansion rate is set by the matter density and once the local curvature 
falls below Rq the expansion rate gets effect of gravity modification. 

For pressure less dust, the effective potential has an extremum at 

2/ - Rf = SirGp. (14) 

For a viable late time cosmology, the field should be evolving near the minimum of the effective potential. The finite 
time singularity inherent in the class of models under consideration severely constraints dynamics of the field. 
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FIG. 1: Plot of the scalaron potential V versus (j) for n — 2 and A = 1.2. The red spot marks the initial condition for evolution 
of <ji. 



The finite time singularity and fine tuning of parameters 



R[t] 




FIG. 2: The evolution of R{t) in the case n = 2 and A = 1.2. The initial condition is as marked by the red spot in figure [T] 

The effective potential has minimum which depends upon n and A. For generic values of the parameters, the minimum 
of the potential is close to (j) — 0, corresponding to infinitely large curvature. Thus, while the field is evolving towards 
minimum, it can easily oscillate to a singular point. However, depending upon the values of parameters, we can 
choose a finite range of initial conditions for which scalar field (j) can evolve to the minimum of the potential without 
hitting the singularity. For instance, let us begin with the initial conditions as shown in figure [T] given by = —0.5 
{R = 1.089). We choose H = 1 then H becomes —1.818 for a given value of R (see eq. (US]))- The minimum of the 
potential for n = 2 and A = 1.2 corresponds to Rmin = 2.218 or 0mm = —0.051 which is indeed close to singularity. 
The numerical results are shown in figures [1] and [21 



The problem of large n, A 



We find that the range of initial conditions allowed for the evolution of (j) to the minimum without hitting singularity 
shrinks as the numerical values of parameters n and A increase. This is related to the fact that for larger values of n 
and A, the minimum fast moves towards = 0, see figure [3l 

For example, in the case of n = 1, A = 4, we find that (j)min = —0.016 which is close to zero. The initial value of (j) 
which is quite near to minimimr also gives the divergence of H{t) and R{t) at the finite time (see Fig. |4|). 
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FIG. 3: Plot of (pmin versus A for different values of n. With increase in n, (pmin moves towards zero (corresponding to infinitely 
large value of R) for smaller values of A. The curves from bottom to top correspond to n = 1, 2, 3, 4 respectively. 

On the other hand, one can see from the figure [3] that also shifts to zero when n is increased. For n — 4 and 

A = 1.1 we have (j)min = —0.003. The range of initial conditions for which the scalar field can evolve to the minimum 
is very small. 



Adding aR^ term and curing the singularity 



We know that in case of large curvature, the quantum effects become important leading to higher curvature corrections. 
Keeping this in mind, let us consider the modification of Starobinsky's model. 



fiR) = R+—R^ + RaX 

Uq 



-1 



1 



then (t> becomes 



2a 



2nX 



(1 + f )"+^ 



(15) 



(16) 



When \R\ is large the first term which comes from aR^ dominates. In this case, the curvature singularity R = ±oo 
corresponds to (j) — ±oo. Hence, in this modification, the minimmn of the effective potential is separated from the 
curvature singularity by the infinite distance in the 0, V{(j)) plane. 



For n — 2, (j) and V{(j)) are given by 



4>{y) = 2ay^ 

V 
Ro 



4Ay 



(l + y2)3 
1 



480(1 + ?;2)6 

- 3(1^^^2)3 (1 + + " (3 + + Ay^)] + \ay^ 

+ ^(32a-7A) tan-i(y). 



(17) 



{X^y (-105 - 595y2 - 2154?/ + 106?;^ -I- 595y** + lOSy^")} 



(18) 



For n = 2, A = 2 and a = 0.5, we have a large range of the initial condition for which the scalar field evolves to the 
minimum of the potential. Though the introduction of R^ term formally allows to avoid the singularity but can not 
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alleviate the fine tuning problem as the minimum should be brought near to the origin to respect the solar constraints. 
Last but not the least one could go beyond the approximation (see equation (HU) by iterating the trace equation 
and computing the corrections to R given by equation (fT4| . As pointed by Starobinsky[i^, such a correction might 
become large in the past. This may spoil the thermal history and thus needs to be fine tuned. The aforesaid discussion 
makes it clear that HSS models are indeed fine tuned and hence very delicate. 

R[t] 
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FIG. 4: The evolution of R{t) forn = 1 and A = 4. 




FIG. 5: Plot of the effective potential for n — 2, X = 2 and a = 1/2 in presence of correction. The minimum of the effective 
potential in this case is located at </>min = 3.952 {Rmin = 3.958). 



III. PARAMETRIZATION AND ANALYSIS 



The Friedmann equation can be casted in the form, 



^2 ^ &^P+Rf'-f 



dif' + RNf") 



(19) 
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where TV = In a. To parameterise the system in a more convenient form we define: 



R 



Since \R\ ^ Rq we have a; ^ 1. The /, /' and /" for the Starobinsky model are of the form 

/ 



7^2 



/' 



X 



2nX 



X 

-(l+n) 



X- 

(2n+l) 







1(- 







(20) 



(21) 



The trace equation can be re-written as 



- + 2A 

a; 



(22) 



Here ftm — f2mo(l + z)^ is the matter density of the universe at a given redshift. f2„iOi is present fractional matter 
density and flA, is the present fractional density of the vacuum. 



The Friedmann equation (equation (|19p) can now be written as 



^0^ 



{xf - SflAH^f") 



(23) 
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Where f{R) and its derivatives are as given above. Note that we can recover the usual Friedmann equation for 
f{R) = R, the case for which our action reduces to Einstein Hilbert action. 

In this paper we test the viability of f{R) cosmology. The Friedmann equation evaluated at z = imposes a constraint 
on the free parameters n and A. 

To check the compatibility of Starobinsky model with observations we need to first obtain the Hubble parameter H, 
as a function of redshift z. The Hubble equation in its present form depends on the curvature i?, through / and its 
derivatives. Therefore, we need to solve the trace equation (equation ([22| ). relating curvature with redshift. In order 
to study the behavior of trace equation, which does not admit a simple solution for general n, we first analyse its 
quantitative behavior for large values of z, corresponding to a; <C 1. Therefore, the trace equation can be simplified 
considerably and in the leading order we recover 



(24) 



(25) 



Similarly, the Friedmann equation in the leading order gives 

Thus at large redshift we recover the standard ACDM cosmology for A = 2, independently of n (see figure [S]). 

There is no general solution for the trace equation and therefore we analyse this model for fixed values of n. It was 
shown by Capozziello and Tsujikawa [l^l that we need n > 0.9 to satisfy solar system tests. In this work we consider 
model with n = 1 and 2. 
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A. n=l 

The trace equation in this case simplifies to a quintic equation 



F{x) 



A 4 
X 



2x^ 



2^x- 

li ^777, 



x\l 



2X^x'\ 



= 0, 



(26) 



where fi^ — $7^0(1 + z)'^- This equation cannot be solved analytically for general x in terms of other parameters. 
Since x <C 1, we approximate this equation with a quartic equation. This approximation is found to be sufficiently 
accurate for a; ^ 1 (see Fig. 




FIG. 6: Plot of F(x) 
dotted line). 



at 2: = 0, A = 2 (indicated by solid line) vs truncated expression of F(x) at 2: = 0, A = 2 (indicated by 



We need to solve for x (x = Rq/R) using the trace equation. We have two quartic equations F{x, z,X) =0 and 
G{x,z,X) = for n =1 and 2 case, respectively. We shall solve the quatric equation with ^Imo = 0.3 and f^A = 0.7. 
Let us first consider the case of n = 1 and look for the roots of = at z = 0. We numerically analyse the Eq. 26 
and find that for generic values of A (0 < A < 10), two of the four roots of quatric equation are always imaginary. The 
Friedman equation puts a constraints on the values of A such that H{z) Hq for z — > 0, see Eq. 23. This condition 
fixes the value of A to 2 for one of the roots, the other root corresponds to A ~ 2.44. We should now check the viability 
of the roots by invoking the non-zero values f the redshift. As z increases, i.e., we move to past, x = Rq/R should 
decrease as R should increase during a viable cosmic evolution which happens in case of A = 2. However, for the root 
corresponding to A ~ 2.44, x increases in the past and x quickly goes beyond its normal range < x < 1 and we 
need not consider this root at all. Therefore we have only one root obtained at A = 2 which produces the correct final 
asymptotic state as a de Sitter model when z goes to -1 (see figure [7]). In this figure we plot the effective equation of 
state with redshift: 



Weff 



2(1 + z) dH 
3H d^' 



(27) 



Once A is fixed there is no other free parameter in the theory. The Hubble parameter can now be plotted for the best 
fit value of A. 
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FIG. 7: Effective Equation of State (EOS) vs. 1+z. Solid line corresponds to n=l and the knotted line corresponds to n=2. 



B. n=2 

The trace equation in this case is a seventh degree polynomial in x 

G{x) ^x'~ + _ 3 ^4 ^ 3 A ^ 2A^) - 3^x' + 2X^) ^ - ^ ^ 0. (28) 

Again we can approximate this equation by a quartic equation in x using the fact that x ^ 1 (see Figl5]). In order 
to calculate the roots of this equation, we followed the same procedure as described for n = 1 case. Here also A = 2 
reproduces the correct late time behavior (see figure [7]). 



IV. OBSERVATIONAL CONSTRAINTS 



In this paper, we check the compatibility of the model with the H{z) data and the measurements of baryon acoustic 
oscillation peak. We find the values of to see how well the /(i?) model with n — 1 and n = 2 accommodate the 
observations. Here, for a given model, we calculate using H(z) data, BAO measurement and for the joint data 
(H(z) + BAO). 



Constraints from H(z) data 

Simon, Verde and Jimenez (2005) employed differential ages of passively evolving galaxies to get the Hubble parameter 
as a function of redshift, H{z) |25|] . The nine data points of H(z) with 0.09 < z < 1.75 have been obtained by using 
absolute ages of 32 galaxies taken from the Gemini Deep Deep Survey (GDDS) and the archival data. 

For the adopted f{R) model, we calculate the values of for n = 1 and n = 2. For this we define 



i=l 



(29) 
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FIG. 8: Plot of G(x) at z = 0, \ — 2 (indicated by solid line) vs truncated expression of G(x) at z = 0, A = 2 (indicated by 
dotted line). 



Here Hq is the present day value of the Hubble constant, Hcxp{zi, n, Hq) is the expected value of the Hubble constant 
in the /(i?) cosmology at redshift Zi for a particular n and Hq. iJobs is the observed value and is the corresponding 
1 a uncertainty in the measurement. The sum is over all observed data points (nine in number). 

As the value of is highly sensitive to the value of Hq, we marginalize over Hq to obtain the modified x^- For this 
we define the likelihood function as: 

Here P{Ho) is the prior probability function for Hq which is Gaussian: 



P{Ho) - exp 



1 {Hp - ggb^)^ 
2 



'Ho 



with Hq^^ as the value of Hq (and is the error in it) as suggested by independent observations. In this paper, we 
also study the effect of different priors on the result. We use two set of priors: 

Set A: Hl^^^ = 68 ± 4 Km/s/Mpc, as obtained from the median statistics analysis of 461 measurements of Hq [2^ . 
Set B: Hq^"" = 77± 4 Km/s/Mpc, as suggested by the Chandra X - ray Observatory results 

We observe that the value of the modified Xh(z) sensitive to the choice of the prior. 

Baryon Acoustic Oscillation (BAO) 

Before recombination, the universe was in the completely ionized state. The cosmological perturbations in the rel- 
ativistic primordial baryon-photon plasma produced acoustic oscillations. These oscillations were imprinted in the 
form of peaks in the late time power spectrum of non-relativistic matter. This acoustic peak which is predicted at 
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FIG. 9: Variation of H{z) with z, for Ho = 77 Kms ^Mpc ^ . Upper curve corresponds to n=l case and the lower curve is for 
n=2 model. The points (boxes) on the lower curve correspond to ACDM. 

the measured scale of 100 Mpc is detected in the large scale correlation function of 46478 sample of luminous red 
galaxies in the Sloan Digital Sky Survey [H]. It is described by a dimensionless parameter A, 



A{n) 



n 



1/2 



r2(z*,n) 



1/3 



(30) 



where = 0.35, r(z,, n) = J^" dz/£{z^, n) is the dimensionless comoving distance to z*, £(z*, n) is given by H{z)/Ho. 

We define Xbao — (-^l"-) ~ ^obs)^/o'3i, with Aohs = 0.469 ± 0.017. The values of Xbao for n = 1 and n = 2 are listed 
in Table 1. 



n Ho prior Xh{z) Xbao Xxotai 

1 68±4 7.29 0.54 7.83 0.87 

1 77±4 6.63 0.54 7.17 0.8 

2 68 ±4 7.79 0.88 8.67 0.96 
2 77 ±4 6.30 0.88 7.18 0.8 



TABLE I: Here Xh 



(^) 



i ln(£). 



V. RESULTS AND DISCUSSION 



The observed late time acceleration of the universe is one of the major unsolved problem in the cosmology. It may 
hint at the breakdown of Einstein GR . This has lead to modification of the Einstein theory of gravity. One of the 
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attractive possibility to modify this theory is to replace the Ricci scalar R with the generic function f{R) in the 
Hilbert action. 



f{R) theories are usually studied by two methods: metric and Palatini approach. The Palatini approach is explored 
extensively in the literature both theoretically and observationally [2^( see the other ref's given in these papers). This 
formulation gives second order differential field equation which can explain the late time behavior of the universe. 



The metric approach of /(i?) theory leads to a fourth order non linear differential equation in terms of scale factor. 
This equation is difhcult to solve both analytically and numerically even for the special cases. Therefore, not much 
observational tests have been performed on f{R) theories based on metric formulation. Our work is an attempt to 
check the concordance of f{R) theory of gravity using the metric approach with some of the cosmological observations. 
In particular, we have explored the Starobinsky model in which function f{R) is analytic, satisfying the condition 
/(O) = 0. This model also passes the Solar system and laboratory tests successfully for large values of n. In this work 
we study the f{R) model for n = 1 and n = 2. 



In order to study the cosmological viability of this theory, we investigate this class of model with two observational 
tests. The first method is based on the Hubble parameter versus redshift data, H{z). The Hubble parameter is related 
to the differential age of the universe through this form 



H{z) 



1 dz 
1 + z dt' 



By estimating the dt/dz, one can obtain directly the Hubble parameter, H{z) at different redshifts. The H(z) data has 
one major advantage, unlike in the standard candle approach (SNe la): the Hubble function is not integrated over. 
The other important feature of this test is that differential ages are less sensitive to systematic errors as compared 
to the absolute ages [1^. This observational H{z) have been used earlier also to constrain various other dark energy 
models 1301. 



In this work, we check the compatibility of f{R) model with n — 1 and n — 2 with the H{z) data. Since Hq is a 
nuisance parameter, we marginalise over Hq. We further combine the results obtained from H{z) data set with the 
BAO data. To perform the joint test, we define the quantity: 

Xtotal — Xh(z) + XbAO- 

(31) 

The results are given in the Table 1. It appears that n = 2 is favored by the observations. This fact is also in 
agreement with solar and laboratory test. The other important conclusions are as follow: 

• The variation of H{z) with redshift becomes independent of n after the redshift around z — 1.8, see Fig. 7. 
Before this redshift ( z < 1.8), there is a small difference between the behavior of H{z) with z for n = 1 and 
n = 2 models. This variation is still within the error bars of the data points. At higher redshift ( i.c when 
R >> Rq) we recover the standard ACDM universe for A = 2, for all the values of n. Hence the thermal history 
of the universe is correctly reproduced by this model. 

• The expansion history for this model with n = 2, A = 2 matches exactly with the standard cosmological model 
(ACDM) with equation of state parameter ui — — l,r2m ~ 0.3 and Ha — 0.7 (see fig. 7). The variation in 
parameters flm and Qa (keeping flmo + S^a — 1) leads to very small changes in values of A, keeping it close to 2. 

• As shown in Table 1 , it seems that the present observational data used in this work prefer n = 2 over n = 1 , as 
indicated by per degree of freedom, xt- However, n = 1 cannot be ruled out. 

• As stated earlier, the algebraic expressions become cumbersome for larger values of n. It is quite possible that 
observational tests discussed in this paper may be compatible with values of n larger than two which are also 
consistent with solar and laboratory tests. However, it should be emphasized that for a given value of A, the 
minimum of the scalaron potential gets closer to zero for larger values of n making the model vulnerable to 
singularity, see Fig. 3. 



Due to large error bars in the data sets used, the variation of the expansion history of the universe studied in this 
paper can easily be accommodated by the observations, thereby making the models close to ACDM at the background 
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level. At present the sample of H{z) data is too small and the error bars are large. In the future, large amount of 
precise H{z) data is expected to become available. This will not only reveal the fine features of the expansion history 
of the universe but also tightly constrain the cosmological parameters. 

There is a further need to explore this model with time based observational tests, like age of the universe and high 
redshift objects. It is known that the evolution of the age of the universe with the redshift vary from model to model. 
So it is possible that the model of the universe which are able to reproduce the total age of the universe at z = 0, may 
not accommodate the objects at high redshift [31| . Therefore the time based observational test may play the key role 
in short-listing the viable dark energy model in the near future. It is also expected that the study of matter power 
perturbations would allow to distinguish this model from the standard cosmology. 
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